Mott scattering at the interface between a metal and a topological insulator 
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We compute the spin-active scattering matrix and the local spectrum at the interface between 
a metal and a three-dimensional topological band insulator. We show that there exists a critical 
incident angle at which complete (100%) spin flip reflection occurs and the spin rotation angle 
jumps by n. We discuss the origin of this phenomena, and systematically study the dependence of 
spin-flip and spin-conserving scattering amplitudes on the interface transparency and metal Fermi 
surface parameters. The interface spectrum contains a well-defined Dirac cone in the tunneling 
limit, and smoothly evolves into a continuum of metal induced gap states for good contacts. We 
also investigate the complex band structure of Bi2Se3. 

PACS numbers: 73.20.-r, 75.70.Tj, 85.75.-d 
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I. INTRODUCTION 

Recently discovered three dimensional topological 
band insulators 1-3 , such as Bii-^Sb^ 4 and Bi2Se3 5 ~ 7 , are 
spin-orbit coupled crystal solids with a bulk gap but pro- 
tected gapless surface states. The low energy excitations 
at the surface are helical Dirac fcrmions, i.e., their spin 
and momentum are entangled (locked) 8 . The charge and 
spin transport on the surface of a topological insulator 
are intrinsically coupled 9 . This makes these materials 
a promising new platform for spintronics. In addition, 
heterostructures involving topological insulator, super- 
conductor, and/or ferromagnet have been predicted to 
show a remarkable array of novel spectral and transport 
properties (for review see Ref. 10 ~ 12 ). 

Electronic or spintronic devices based on topological 
insulators will almost inevitably involve metal as mea- 
surement probes or functioning components 13 . This mo- 
tivates us to study the local spectrum near the interface 
between a metal (M) and a topological insulator (TI). 
For a metal-ordinary semiconductor junction with good 
contact, it is well known that the metallic Bloch states 
penetrate into the semiconductor as evanescent waves 
localized at the interface (for energies within the band 
gap). Such interface states are known as metal induced 
gap states (MIGS) 14 ' 15 . They play an important role 
in controlling the junction properties, e.g., by pinning 
the semiconductor Fermi level to determine the Schottky 
barrier height 16 , a key parameter of the junction. 

The local spectrum at the M-TI junction is intimately 
related to the spin-active scattering of electrons at the M- 
TI interface. In this paper, we systematically study the 
evolution of the scattering matrix and the interface spec- 
tra with the junction transparency and metal Fermi sur- 
face parameters. The scattering matrix 17 we obtain here 
also forms the basis to investigate the details of the super- 
conducting proximity effect near the superconductor-TI 
interface 18 , which was shown by Fu and Kane to host 
Majorana fcrmions 19 . 

The scattering at the M-TI interface differs signifi- 



cantly from its two dimensional analog, the interface be- 
tween a metal and a quantum spin Hall (QSH) insulator 
studied by Tokoyama et al 13 . They predicted a giant 
spin rotation angle a ~ ir and interpreted the enhance- 
ment as resonance with the one-dimensional helical edge 
modes. By contrast, for M-TI interface we predict a crit- 
ical incident angle at which complete spin flipping occurs 
and the spin rotation angle jumps by n. We will explain 
its origin, in particular its relation to the surface helical 
Dirac spectrum, and discuss its spintronic implications. 
This paper is organized as follows. We will first com- 
pute the scattering matrix using a k • p continuum model 
by matching the envelope wave functions at the M-TI 
interface. This simple calculation is easy to understand, 
and it brings out the main physics of our problem. Along 
the way, we will discuss the complex band structure of 
Bi2Se3, which describes the decaying (rather than propa- 
gating Bloch wave) solutions of the crystal Hamiltonian. 
The various caveats of this calculation are then reme- 
died by considering a much more general lattice model. 
Most importantly, it enables us to track how the scatter- 
ing matrix and interface spectrum change with interface 
transparency. It also sheds light on the origin of perfect 
spin-flip scattering at the critical angle. We will show 
that the results obtained from these two complementary 
methods are consistent with each. 



II. MODEL HAMILTONIAN AND COMPLEX 
BAND STRUCTURE 

We consider Bi 2 Se3 as a prime example of 3D strong 
topological insulators. Its low energy k • p Hamiltonian 
was obtained by Zhang et al 6 , 
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FIG. 1: (a) Scattering geometry at a metal (M)-topological 
insulator (TI) interface, (b) Schematic band structure of the 
metal (modeled by Hm) and topological insulator. 
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FIG. 2: The complex band structure of topological insulator 
described by H T i(k) for k y = 0, k x = 0.02 (left) and 0.04 
(right). E is measured in eV, and k in A -1 . Subgap states 
with complex k z represent evanescent waves. The topology of 
real lines 25 changes as k x is increased. 
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2(kx + ky). The numerical values of M, A, B, C , D are 
Ref. 6 . 



given in Ket". We choose the basis (|+ f), |+ I), | — t 
),|— I)), where ± labels the hybridized p z orbital with 
even (odd) parity 6 . The Gamma matrices are defined 
as To = T3 ® 1, Ti = fi ® &i, with fj (<7j) being the 
Pauli matrices in the orbital (spin) space. The chemical 
potential of as-grown Bi 2 Se3 crystal actually lies in the 
conduction band 8 . By hole doping 8 or applying a gate 
voltage 20 , the chemical potential can be tuned inside the 
gap. The system is well described by Hti (note that 
energy zero is set as in the middle of the band gap) . 

In this section, we first adopt a rather artificial model 
for metals with negligible spin-orbit coupling. It is ob- 
tained by turning off the spin-orbit interaction (setting 
dn = for /i=l,2,3) in Hti and shifting the Fermi level 
into the conduction band. The result is spin-degenerate 
two-band Hamiltonian 

H M {k) = [e (k) -E F }1 + do(k)f 0. 

Its band structure, schematically shown in Fig. 1(b), 
consists of two oppositely dispersing bands (the solid and 
dash line). Ep is tuned to be much higher than the band 
crossing point, so the scattering properties of low energy 
electrons near the Fermi surface are insensitive to the 
band crossing at high energies. This claim will be verified 
later using a more generic model for the metal. A similar 
model was used in the study of metal-QSH interface 13 . 

Matching the wave functions of two dissimilar mate- 
rials (such as Au and Bi 2 Se3) at interface is in general 
complicated within the k • p formalism, because the en- 
velope wave functions on either side are defined using 
different basis (see Ref. 21 and reference therein). For 
the particular model Hm , however, such complication is 
circumvented. Then, then wave functions at the metal- 
TI interface (z — 0) satisfy the Ben-Daniel and Duke 
boundary condition 22 , 

$M = &TI, i>M$M = VTI&TI- 

Here <&i is the four-component wave function, and the ve- 
locity matrix Vi = dHi/dk z , i e {M,TI}. Such bound- 



ary condition assumes good atomic contact between two 
materials. 

We are interested in energies below the band gap of TI, 
so <$>ti is evanescent in nature and only penetrates into 
TI for a finite length. Such localized (surface or interface) 
states inside topological insulator can be treated within 
the k • p formalism using the theory of complex band 
structures, pioneered by Kohn 23 , Blount 24 , and Heine 25 
et al. The main idea is to allow the crystal momentum 
to be complex and analytically continue iJjv(k) to the 
complex k plane. While the extended Bloch waves are 
the eigen states of ifr/(k) for real k, eigen functions of 
Hti{^) for complex k describe localized states. Together 
they form a complete basis to describe crystals of finite 
dimension. 

In our scattering problem, we have to find all eigen 
states of J?T/(k) with energy E and wave vector k = 
(k x , k y , k z ), where k x and k y are given and real, but k z is 
complex and unknown. For a general k • p Hamiltonian 
such as Hti , we follow Chang and Schulman 26 to rewrite 
it as 

Hti = h (k x ,k y ) + h x k z + h 2 k 2 z , 

where hi = A1D3, and h 2 = —BiTq. Then the eigen 
equation (Hti ~ El)(f> = can be reorganized into an 
eigen value problem for k z , 
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Then all possible values of k z can be obtained for given 
incident parameter E, k x , and k y . For the anisotropic 
Dirac Hamiltonian HtiQ<.), the energy eigenvalues can 
be obtained analytically 28 , which allows for an analytical 
solution of the complex band structure. 

For E within the gap, there are in general 4 pairs of 
complex solution of k z , for if k z is a solution so is k*. 
We label those with positive imaginary parts with {k z }, 
and the corresponding wave function {4> u }, v = 1, 2, 3, 4. 
They are decaying solutions in the half space z > 0. 



In our model, k z turns out to be doubly degenerate, as 
shown in Fig. 2. The wave function inside TI (z > 0) 
then has the form 



*TI = $>" e< ** 



III. SCATTERING MATRIX FROM 
WAVE-FUNCTION MATCHING 

To set the stage for discussing scattering off a topolog- 
ical insulator, it is instructive to recall the generic fea- 
tures of elastic scattering of electrons by a heavy ion with 
spin-orbit interaction. This classical problem was solved 
by Mott, and known as Mott scattering. The scattering 
matrix has the general form 27 



''Mott 



ul 



(k-i x k j. 



where kj and k are the incident and outgoing momen- 
tum respectively, & is the Pauli matrix, and u, w de- 
pend on the scattering angle. It is customary to de- 
fine the spin- flip amplitude / = 5 2 i, and spin-conserving 
amplitude g = Si±. Both / and g are complex num- 
bers, their relative phase defines the spin rotation angle 
a = Arg(<7*/). One immediately sees that for back scat- 
tering, Smou = ul, so there is no spin flip, / = 0. As we 
will show below, this also holds true for scattering off TI. 
Now consider an electron coming from the metal with 
momentum k incident on the M-TI interface located at 
z = 0, as schematically shown in Fig. 1(a). We assume 
the interface is translationally invariant, so the transverse 
momentum ky = (k x , k y ) is conserved, and the energy E 
of the electron lies within the band gap of TI. Then, only 
total reflection is possible, but the spin-orbit coupling in- 
side TI acting like a k-dependent magnetic field rotates 
the spin of the incident particle. The scattering (reflec- 
tion) matrix has the form 



S(k) 



9 I 
f 9 



where |ff| 2 + |/| 2 = 1- Our goal is to find the dependence 
of the scattering amplitudes /, g on k, or equivalcntly, 
on energy E and incident angle 9. From time-reversal 
symmetry, f(E,6) = f(E,-0) and g(E,6) = g(E,-9). 
We shall show that /(ky) = -/(-ky), fl (k,|) = ,g(-k||). 
So / is an odd function of 6, while g is even in 6. Since our 
problem can be viewed as coherent multiple scattering 
from a lattice array of Mott scatters occupying half the 
space, we will refer to spin-active scattering at the metal- 
TI interface as Mott scattering. 

Consider a spin up electron from the conduction band 
of the metal with momentum k and energy E = eo(k) — 
Ep — rfo(k) lying within the band gap of TI. The wave 
function inside the metal (z < 0) has the form 



up to the trivial e l ^ kxX+kyV ^ and renormalization factor. 
Here k z = z ■ k, and {r^} are the reflection amplitudes. 
We identify the spin flip amplitude / = r^ and the spin- 
conserving amplitude g = r$. Note that there is no prop- 
agating mode at energy E available in the valence band 
for the reflected electron. So k' z is purely imaginary. At 
such energy, there is no propagating mode available in 
TI. We have discussed the evanescent wave function $jv 
in the previous section. With $m and 4>ti , we solve the 
boundary condition at z = to obtain r Vl t„ and the 
scattering matrix S. 

Fig. 3 shows the magnitude and phase of / and g ver- 
sus the incident angle 9 for E — O.leV, with Ep set to be 
0.28eV. At normal incidence, 6 — 0, spin flip scattering 
is forbidden as in the single-ion Mott scattering. With 
increasing the magnitude of g drops continuously. At 
a critical angle 9 C1 \g\ drops to zero and we have perfect 
(100%) spin flip reflection. At the same time, the spin 
rotation angle a (the relative phase between / and g) 
jumps by tt. 

It is tantalizing to think of what happens at 9 C as res- 
onant scattering with the helical surface mode of the TI. 
This however is problematic. We are considering good 
contacts at which the wave functions of the two mate- 
rials hybridize strongly. Surface mode is preempted by 
MIGS. Indeed, we checked that the corresponding crit- 
ical transverse momentum ku depends only weakly on 
E. This is at odds with the linear dispersion of the TI 
surface mode, E = A 2 k^ . To gain better understand- 
ing, we now switch to a lattice model to systematically 
study the role of interface transparency and metal Fermi 
surface parameter (Ef,kf,v/) on the scattering matrix. 



IV. INTERFACE SPECTRUM AND 

SCATTERING MATRIX FROM LATTICE 

GREEN FUNCTION 

We consider a simple lattice model for the M-TI junc- 
tion. The topological insulator is modeled by a tight 
binding Hamiltonian on cubic lattice, 



,71+1 



h.c. 
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M 



-ik^z Ak z z 



(r ie -^ z , r 2 e~ lk * z , e lk * z + r 3 e-^ z , r 4 e- ik * z ) T , 



k||,n 

+ ^k|,,n M( k l|)ro +a 2 (t 1 sin k x +t 2 smk y ) tpk^nt ■ 

Here ip = (ip+^,?p+i,ip-^, ip-i) T is the annihilation op- 
erator, d(k||) — M — 2b i + 2b 2 (cosk x + cosk y — 2) with 
k measured in 1/a. The cubic lattice consists of layers 
of square lattice stacked in the z direction, n is the layer 
index, and ky is the momentum in the xy plane. The 
isotropic version of Jtf'pi, with a\ — a 2l b\ — b 2 , was stud- 
ied by Qi et al as a minimal model for 3D topological 
insulators 28 . To mimic Bi 2 Se3, we set the lattice spacing 
a = 5.2A, which gives the correct unit cell volume, and 
otj = Aij a, bi = Bi/a 2 for i = 1, 2. Although a crude car- 
icature of the real material, M'pi yields the correct gap 
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FIG. 3: The magnitudes (upper panel) and the phases (lower 
panel) of the spin-flip amplitude / and spin-conserving ampli- 
tude g versus the incident angle 0. E = O.leV, i?F=0.28eV. 
\d\ 2 + l/| 2 = 1- Arg(<?) and Arg(/) are shifted upward by n 
for clarity. 
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FIG. 4: The spectral function N(E,k x ,k y — 0) at the inter- 
face of metal and topological insulator. Left: good contact, 
J = £m, showing the continuum of metal induced gap states. 
Right: poor contact with low transparency, J = 0.2tM, 
showing well defined Dirac spectrum as on the TI surface. 
tu = 0.18eV, Hm = — 4tjv/, a is lattice spacing. 



size and surface dispersion, it also reduces to the contin- 
uum k • p Hamiltonian Hti in the small k limit, aside 
from the topologically trivial eo(k) term. 

As a generic model for metal, we consider a single band 
tight binding Hamiltonian on cubic lattice, 



"/S ! 






[/i(k|i)n k|l 



t M 4>L 



h.c. 



where /i(ky) = — 2tM(cosk x + cosk y ) — hm- The Fermi 
surface parameters of the metal can be varied by tuning 
£m and \im- The metal occupies the left half space, n < 
0, and the TI occupies the right half space n > 1. The 
interface domain consists of layer n = 0, 1. The coupling 
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FIG. 5: The spin-conserving reflection amplitude \g\ and 
spin rotation angle a versus the incident angle 9 for increas- 
ing contact transparency, J jtu = 0.25, 1, 1.5,2 (from left to 
right). t M = 0.18ey, fj, M = -4i M , E = 0.05eV, k y = 0. 
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between metal and TI is described by hopping, 

^LR = ~ Yl J ^k| h n=l/,<T^k || ,n=0, CT +^-C. 
k|| ,£,cr 

Ji is the overlap integral between the p-orbital I — ± 
of TI and the s-like orbital of metal. For simplicity, we 
assume Ji is independent of spin. Then, J + = — J_ = J. 
J can be tuned from weak to strong. Small J mimics a 
large tunneling barrier between M and TI, and large J 
(comparable to tu or B-i) describes a good contact. 

The lattice Green function of the composite system 
is computed via standard procedure by introducing the 
inter-layer transfer matrix and the method of interface 
Green function matching 29 . Fig. 4 shows two examples of 
the local spectral function (momentum-resolved density 
of states) at the interface, 



7V(E,k|, 



^ imTrSf(£,k||) r 



n=0,l 



where < ^(E,'k^) nin > is the local Green function at the in- 
terface with n, n! = 0, 1, and the trace is over the spin and 
orbital space. In the tunneling (weak coupling, small J) 
limit, the interface spectrum includes a sharply defined 
Dirac cone as on the surface of TI. As J is increased, the 
linearly dispersing mode becomes ill defined and eventu- 
ally replaced by a continuum of metal induced gap states. 
Once the lattice Green function is known for given in- 
cident E and fcn , the scattering (reflection) matrix can 

be constructed from Sf by 29 , 

S(£,k||)=^(£,k,|) ,oS M ^k||)-i 

where gu is the spin-degenerate bulk Green function of 
metal. Fig. 5 shows the evolution of \g{6)\ and a(0) 
for increasing J, where a level broadening of E/10 is 
used. Most importantly, we observe that the existence of 
a critical angel 8 C , where complete spin- flip occurs and a 
jumps by 7r, is a robust phenomenon. It is independent 
of the details of the contact, the metal Fermi surface, or 
other high energy features in the band structure. 



To understand the perfect spin flip, we first focus on 
the tunneling limit, J <C £m- In this limit, the local 
spectrum at layer n = 1 as shown in the right panel of 
Fig. 4 approaches the TI surface spectrum, namely the 
helical Dirac cone. An incident up spin tunneling across 
the barrier will develop resonance with the helical mode, 
which is a quasi-stationary state with long life time, if its 
momentum and energy satisfy k» = E/A^. Moreover, it 
has to flip its spin, since only down spin can propagate 
in the k x direction (suppose k y — 0). The tt jump in the 
phase shift is also characteristic of the resonance. Indeed, 
we have checked that precisely at 9 C the resonance cri- 
terion, kf sin9 c — E/A 2 , is met. We also varied \xu for 
fixed J and £m, bigger \iu yields a bigger Fermi surface 
and a smaller 9 C . This is consistent with the resonance 
criterion above. 

As J is increased, the width of the resonance grows 
and eventually it is replaced by a broad peak (dip) in 
|/| (|g|), but the vanishing of \g\ and ir shift in a at 9 C 
persist to good contacts, even though in this limit the 
interface is flooded by MIGS (left panel of Fig. 4) and 
bears little resemblance to the Dirac spectrum. With all 
other parameters held fixed, 9 C increases with J. Quali- 
tatively, coupling to TI renormalizes the metal spectrum 
near the interface, producing a smaller effective kf (hence 
a larger 9 C ) compared to its bulk value. It is remarkable 
that perfect spin flip at the critical angle persists all the 
way from poor to good contacts. Indeed, the main fea- 
tures observed here for for good contacts using the lattice 
model agree well with the results obtained in previous 
section by wave function matching. 



by ARPES (or scanning tunneling microscope) experi- 
ments on metal film coated on a topological insulator. 
Our results also suggest that a topological insulator can 
serve as a perfect mirror to flip the electron spin in metal. 
Such spin-active scattering at the M-TI interface may be 
exploited to make novel spintronic devices. The magni- 
tude of g or / can be measured by attaching two fer- 
romagnetic leads to a piece of metal in contact with 
TI, forming a multi-terminal device. One of the ferro- 
magnetic leads produces spin-polarized electrons incident 
on the M-TI interface at some angle, while the other 
lead detects the polarization of reflected electron, as in 
a giant magneto-resistance junction. The spin rotation 
angle a can be measured indirectly by comparing the 
predicted current-voltage characteristics of M-TI-M or 
Superconducto-TI-Superconductor junctions, which are 
sensitive the phase shift a. It can also be inferred from 
the spin transport in a TI-M-TI sandwich, as discussed 
for QSH insulator in Rcf. 13 . Detailed calculations of the 
transport properties of these structured, using the scat- 
tering matrix obtained here, will be subjects of future 
work. 
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We now discuss the experimental implications of our 
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